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■ Abstract 



The first four terms of an expansion in M^/Mf of the total inclusive cross section 
'■ for Higgs production in gluon fusion are evaluated through next-to-next-to-leading 

I order QCD. A reliable and precise approximation of the full top mass dependence at 

QQ ■ NNLO is derived and compared to the frequently used heavy-top limit. It is found that 

both results agree numerically to better than 0.5% in the Higgs mass range of 100- 
300 GeV. This validates the higher order results for the inclusive Higgs cross section 
and justifies the heavy-top limit as a powerful tool for Higgs phenomenology at the 
LHC and the Tevatron. 



^ . 1 Introduction 



It is expected that the Large Hadron Collider (LHC) will provide insight into the mecha- 
' nism of electro-weak symmetry breaking. Many of the theoretical models, including the 

CN ■ Standard Model (SM), implement it by introducing an as of yet unobserved elementary 

scalar particle, a so-called Higgs boson. One of the tasks of the LHC is therefore to search 
for a Higgs boson and to measure its properties (for reviews, see Refs. [1, 2, 3]). 



0^ , 

, Precise measurements will require a thorough understanding of the Higgs production cross 

I section. In this respect, the gluon fusion process poses a number of problems: the gluon 

parton densities need to be known very precisely; the cross section is of order and thus 
very sensitive to the precise value of this quantity; the radiative corrections and thus the 
rS I estimated uncertainty from higher order effects is unusually large. 



The latter aspect is particularly problematic because the gluon fusion process is loop 
induced, so that the next-to-leading order (NLO) effects already require a two-loop calcu- 
lation. Fortunately, it was found that the NLO ET-factor is extremely well reproduced in 
an effective theory approach, obtained by integrating out the top quark from the theory. 
Details will be discussed below. 

Based on this observation, it is commonly assumed that the NNLO corrections, obtained 
within this effective theory approach, are practically equivalent to the full calculation, at 
least below the top quark threshold. In fact, almost all radiative corrections beyond NLO 
have been treated in this framework up to now. 

In this paper, we study the quality of this approximation by evaluating subleading terms in 
the expansion in l/Mf. After introducing our notation and the basic formulas in Section 2, 
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we describe our calculational methods in Section 3, while the analytic results for the 1/Aft 
expansion are presented in Section 4. A matching of the low- to the high energy region is 
described in Section 5 and the numerical analysis is presented in Section 6. The conclusions 
of our findings are drawn in Section 7. 



2 Higgs production through gluon fusion 

The inclusive hadronic cross section for SM Higgs production in proton-(anti-)proton 
collisions is obtained by convoluting the partonic cross section dap-^n+x for the scattering 
of parton a with parton [3 by the corresponding parton density functions (pa/pi^) (PDFs): 

app>^H+x{s) = ^ / dT£af3{r)aai3^H+x{s = TS) , 
/■I da; 

£a(i{T)= / [4>a/p{x)(pl3/p'{T/x)] , p'^{p,p}- 

J T 

(In the following, we will shorten or drop the subscripts on £, (p, a, and a whenever there 
is no chance of confusion). In this formula, we have suppressed the dependence on the 
factorization scale np for (p and a, as well as the dependence of a and a on Mt and M^. 
The partonic cross section a is evaluated in terms of a perturbation series in a^. Its 
leading order contribution arises from the triangle diagram shown in Fig. 1 (a) (plus the 
one with opposite fermion direction) , where the fermion line can in principle be any of the 
six quarks. However, due to the Yukawa couplings, the by far dominant contribution is 
due to top quarks, while the bottom quark has an effect of roug ly 7% at the LHC\ and 
9% at the Tevatron. All lighter quarks can safely be neglected. 

NLO QCD corrections to the inclusive cross section have been known for a long time. 
They were first evaluated for the top quark induced terms in the so-called "heavy-top 
limit" which will be described in more detail below, and were found to increase the cross 
section by roughly 70% at the LHC w.r.t. the leading order prediction [4, 5]. This was 
confirmed by the general result for arbitrary top and bottom quark mass [6, 7]. 

The large NLO effects clearly asked for the evaluation of the NNLO corrections which, 
considering its success at NLO, are generally assumed to be well approximated in the heavy- 
top limit. They lead to another significant increase of the total cross section [8, 9, 10], so 
that its actual value turns out to be roughly a factor of two above the LO prediction for 
the LHC, and even up to a factor of three at the Tevatron (for the latest compilations, see 
Refs. [11, 12]). Further studies that go beyond the fixed-order NNLO result have not found 
significant effects and have thus corroborated the stability of the perturbative series (see, 
e.g. Refs. [13, 14]). 

The goal of this paper is to go beyond the heavy-top limit in order to test the quality of 
this approximation. In principle, our approach would also allow us to derive an improved 
prediction of the inclusive cross section. However, it will turn out that the heavy-top limit 
works so well that it reproduces our improved result to better than 0.5% accuracy. 

^By "LHC", we refer to pp collisions at 14TeV center-of-mass energy, although the initial energy of 
the LHC will be lower; more detailed phenomenological studies at various energies are deferred to a future 
publication. 
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We write the top quark induced partonic cross section as 



a/3 ; 



(2) 



with 



vr 



256 



vr 



1 + (1 — r) arcsin 



2 1 



(3) 



Here, Gf ~ 1.16637 • 10 ^ GeV ^ is Fermi's constant, and throughout this paper, Mt 
denotes the on-shell top quark mass, and as = cXs (/^r) the strong couphng for five active 
quark flavours at the renormahzation scale /^r. The kinetic terms assume the form 

2 



Aa/3(x) 



X) + 



(4) 



where 



X 



Ml/s. 



(5) 



The A^") stiU depend on M^, and logarithmically on the renormahzation and factoriza- 
tion scales //R and /xp- At NLO, the full Mt dependence is known in numerical form [6] 
(the virtual terms are known analytically [15, 16, 17]). At NNLO, only the heavy-top 
approximation is known, however [8, 9, 10]. It will be discussed in more detail in the next 
section. 

Note that in order to arrive at a consistent NNLO result for the hadronic cross section, one 
needs to evaluate Eq. (1) not only with an NNLO expression for o", but also by taking into 
account expressions for and the PDFs at the appropriate order. In this paper, we use 
the central set of MSTW2008 [18] which are the latest available NNLO PDFs. A detailed 
study of the PDF uncertainties is left for a future publication. 

Apart from these pure QCD corrections, also the leading electro-weak effects have been 
evaluated [19, 20, 21, 22], and an estimate of the mixed electro- weak/QCD corrections is 
available [12] as well. 



3 Calculation of the 1/Mt terms 

It is well known that the gluon-Higgs interaction in the heavy-top limit can be expressed 
in terms of an effective Lagrangian [23]: 



Cff — L-l <-^/^l/<-^ -I-'Lqcj-), 



(6) 



where v = 246 GeV and the Wilson coefficient Ci is meanwhile known through 0(af ) [24, 
25, 26]. For completeness, we quote it here through 0{a^g) [27, 28] which is sufficient at 
NNLO: 



„ 1 tts f 11 /tts^ ^ 

Ci=--- 1 + --+ - 
3 vr I 4 vr V vr 



2777 19 , 

1 In — ^ + ni 

288 16 M2 



67 1 /i^ 

— In — ^ 

96 3 M? 



(7) 



(5) 

jCqqj-, is the QCD Lagrangian with n; = 5 massless quark flavours. 
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Instead of working strictly within this effective theory, however, one usually factors out 
the full leading order top mass dependence from the Higgs production cross section and 
writes 




where ctq is given in Eq. (3), and Aa/3,00 is evaluated on the basis of Eq. (6). It is then 
clear that does not depend on Alt, while A^-* only has a logarithmic Mj-dependence 
through Ci(as). Note that with this definition of the heavy-top limit, where the full 
Mf-dependence in o"o is kept, the LO cross section is identical to the one in the full theory. 

It was observed long ago [28] that the total cross section at NLO in QCD is approximated by 
cr^^^ to better than 1% up to values of Mh = 2Mf ; even at Mh = 1 TeV, the deviation to 
the exact result remains below 10% (see, e.g., Ref. [29]). This precision is quite remarkable, 
because at higher orders, gluon fusion is a 3-scale process, depending on Mt, Mh, and 
the center-of-mass (cm.) energy While at LO s is fixed to My, the real radiation 
of quarks and gluons starting at NLO allows ^/I to vary from Mh up to the hadronic 
cm. energy which may reach 14 TeV at the LHC. Considering more exclusive quantities 
such as pr-distributions (see Ref. [30] and references therein), or even fully differential 
approaches as in Refs. [31, 32], the number of scales increases further and the heavy-top 
approximation may no longer be under control. Recent studies in this direction can be 
found in Refs. [33, 34]. 

For the inclusive cross section - which is the subject of this paper - it is usually argued 
that the reason for the high quality of the approximation in Eq. (8) is the dominance of 
soft gluon radiation in the higher order corrections, but a solid theoretical justification 
and a quantitative error estimate are still unavailable. One may therefore remain in doubt 
about the applicability of Eq. (8) beyond NLO. A way to test it is to evaluate subleading 
terms in 1/Mt to the partonic cross section, similar to what was done at NLO in Ref. [35]. 
The effective Lagrangian of Eq. (6) could clearly be extended to incorporate such terms 
by including higher dimensional operators (see, e.g., Ref. [36]). However, in general at 
each order in 1/Mt the number of operators grows, and renormalization becomes more 
and more clumsy, for example. 

In this paper, instead of constructing an effective Lagrangian, we directly evaluate the 
Feynman diagrams obtained from the six-fiavour Lagrangian by applying the method of 
asymptotic expansions (for a review, see Ref. [37]). At NNLO, the contributing diagrams 
are at the 3-loop level for the process gg — > H, at the 2-loop level for the single real 
emission processes gg Hg, qg Hq, qq Hg, and at the 1-loop level for the dou- 
ble real emissions gg Hgg, gg — > Hqq, qg — > Hqg, qq — > Hqq, qq — > Hqq (identical 
quark flavours), and qq' Hqq' (different quark flavours; it is understood that the charge 
conjugated processes need to be taken into account as well). Each Feynman diagram con- 
sidered here contains at least one top quark loop. Examples for the various contributions 
are shown in Fig. 1. 

In our approach, we assume Mt heavier than any other basic scale in the process. This 
allows us to express all Feynman integrals as convolutions of massive vacuum integrals with 
at most three loops, and massless 3/4/5-point functions through 2/1/0 loops, respectively. 
In fact, the purely virtual contributions at NNLO have already been calculated using this 
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(d) (e) 




(f) (g) (h) 

Figure 1: (a) LO Feynman diagram for the gluon fusion process; (b)-(h) Sample 
Feynman diagrams contributing to the inclusive NNLO cross section for Higgs 
production in gluon fusion, (b)-(e) single real radiation; (f)-(h) double real radi- 
ation. 
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Figure 2: Diagrammatic representation of the asymptotic expansion of a par- 
ticular Feynman diagram in the hmit s,M^ < AM^. The dia grams left of (S) 
represent subdiagrams of the original diagram that are to be expanded in the 
momenta corresponding to the dotted external lines before the loop integration. 
In this way, it is apparent that the original integral, depending on s, and 
M^, is decomposed into products of "tadpole" integrals with vanishing external 
momenta and massless four-point functions. The shaded blob in the diagrams 
right of (8) represents an effective vertex given by the result of the diagram left of 
(for details of asymptotic expansions, see Ref.[37], for example). 



method [38, 39]^, so we will not discuss them in more detail at this point. 

For the double real emission contributions, the asymptotic expansion is equivalent to the 
interchange of loop integration and Taylor expansion in p/Mt, where p is any component 
of the external momenta, and thus one is left only with 1-loop massive tadpole integrals 
which can be easily evaluated with the help of MAT AD [40] . The difficulty arises from the 
phase space integration which we perform in terms of an expansion around j; = 1 [8] 
("soft expansion"). This is fully justified due to the fact that the 1/Mt-expansion assumes 
< 2Mt and thus x > 0.1 anyway. Apart form that, for the l/Mj^ -terms, it was observed 
that the hadronic cross section converges very well to the full result when successively 
higher order terms in (1 — x) are included in the partonic cross section [8]. 

The graphical representation of the asymptotic expansion of one of the single real emission 
diagrams is shown in Fig. 2. The resulting integrals can be calculated by standard means: 
For the massive tadpoles, we use MATAD [40], and the 1-loop tensor integrals are solved using 
Schwinger parametrisation and integration-by-parts [41]. Both routines are embedded 
in the q2e/exp framework [42, 43] which, in combination with the diagram generator 
qgraf [44], provides a fully automatic way to calculate the relevant diagrams, including 
the procedure of asymptotic expansions. This setup allows for the evaluation of arbitrary 
orders in the 1 /Mt expansion, with the only limitation arising from the available computing 
power. For this publication, we found that the optimal cost/benefit ratio^ is reached at 
0{1/M^), i.e., three orders beyond the heavy-top limit known so far. After multiplication 

■^Compared to Ref. [38], we have additionally calculated the 1/Mf term and found agreement with 
Ref. [39]. 

^As the measure of "benefit" we considered the numerical difference between the highest two terms in 
the 1/Mt expansion. 
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by the corresponding NLO diagrams, we evaluate the phase space integrals in terms of 
hypergeometric functions depending on x = M^/s and e = (4 — D)/2. Expansion around 
e = leads to a Laurent series with poles at e = through 0(l/e^), and coefficients 
depending on x through polylogarithms up to the 4**^ degree. Since we calculated the 
double real emission contributions as expansions in (1 — x), we have to expand the single 
real emission terms in the same way. This is of course trivial from the full x-dependent 
result. However, also here we can perform the soft expansion before the phase space 
integration which serves as a useful check. 

In the sum of the virtual, single-real, and double-real emission contributions, the poles 
cancel, and upon UV-renormalization of as (for which we adopt the MS scheme), Mt and 
the gluon wave function (both in the on-shell scheme), also the 1/e^-terms drop out. The 
remaining poles are of infra-red nature, and are absorbed into the PDFs with the help of 
the usual mass factorization procedure. The relevant convolution integrals are calculated 
by transforming them into Mellin space (which turns them into simple products) and 
subsequent inverse Mellin transformation with the help of Ref . [45] . 



4 Analytical results 



A natural decomposition of the kinetic terms to the inclusive cross section is as follows: 

2n-l 



)<5(l-x)+5]6l")p,( 



fc=0 



(9) 



where the I'fc(x) are the usual plus-distributions defined as 



dxPfc(x)/(x)^ / dx [/(x)-/(l)] 



ln''(l -x) 

1 — X 



(10) 



for arbitrary functions /(x) (differentiable at x = 1). The terms arising from a^"^ and 

(n) 

are called the "soft+virtual contribution", while is referred to as the "hard contri- 
bution" . The a^") , 6^"^ , and h^^l (x) are functions of the top quark mass. We write them 



as 



/ Af H \ * 

{a, hk, hafsix)} = Z^{ -Jf j ^k,i, hai3,i{x)} . 



(11) 



The leading terms in l/Mf, i.e. i = 0, have been obtained in Refs. [46, 47, 48] (soft+virtual) 
and Refs. [8, 9, 10] (hard)."^ 

Concerning the subleading terms in 1/Afj, let us begin with the NLO result (terms with 
odd powers of Mn/Mf vanish): 
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34 



.(1) 



3553 



1- 6C2 , On = , a\ — 

2 ^ 2 135 ^ 113400 



,(1) 



917641 



.(1) 
'^1,0 



^(1) 



190512000 ' (12) 



y i>2, 



''in fact, the leading logarithmic hard term was obtained before in Ref. 
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where C2 = 7r^/6, and we have set /ip = /xr = Mh- For completeness, we have also listed 
the leading terms in l/M^. Note that the coefficients of the plus distributions are fully 
determined in the heavy-top limit, i.e., they do not receive 1/Mt corrections [6, 35]. 6q^q 
becomes non-zero for fip 7^ Mh as is obvious from the evolution equations for fip. 

The hard terms read, for the gg channel: 

(1) _ 3x{l-x) (1) _ (1 - x) (37 + 76x - 83x2 + 68x3) 

V2- ^5 ' KgA- 11200x2 ' 

(i) (1 - x) (3864 - 4008 x + 4133 x^ - 2668 x^ + 3770 x'^) 
h„ 



39'^ 2016000x3 



where 



L(x) = ln(l - x) . (14) 



For the qg channel (which is identical to the qg channel), one finds: 
Co = ^ (2 - 2 X + x2) (2L(x) - In x) - 1 + 2 X - y , 



ll(-4 + 6x-3x2 + x3) 

270 X ' 
24409 - 69264 X + 62052x2 - 24052x3 + 6855 x^ 



(15) 



'i'S'^ 1814400x2 



-181104 + 690137 X - 1008064x2 + 658284x3 - 214280 x^ + 55027 x 



5 



^S'S 108864000x3 
And finally, for the qq channel, it is 

.(1) _ 32 ^3 , (1) _ 88 (l-x)3 _ (1 -x)3 (3487 + 842 x) 

'W.O ~ r,^ V-*- ' ""00.2 ~ ^ntr ^ ' "'a 



h 



'}^'0~27^ ' ' ««"2~ 405 X ' 85050x2 
(1) _ (l-x)3 (41160 + 16271 x + 5573x2) 

9^'^ ~ 5103000x3 



(16) 



Note that except for the terms for the gg channel which can be found in Ref. [35], 

even this NLO expansion has never been given in the literature. For the sake of brevity, 
we quote only the first four terms in here; higher order terms are available from the 

authors upon request. The same is true for the result for general values of //p and //r. 

There are several observations to be made: 

• Pulling out the leading order top mass dependence in terms of o"o, as done in Eq. (2), 
absorbs all the logarithmic x and (1 — x) dependence into the leading l/Mt terms. 

• In general, each power in M^/M^ is accompanied by a power in 1/x. This is a 
consequence of assuming that the top quark mass is the heaviest scale of the process. 
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The expansion therefore is not just in M'^/Mf as suggested by the form of Eq. (11), 
but also in s/Mf which leads to 

A = Mi (17) 

These terms are a sign of the breakdown of the heavy-top limit: for large s, they 
lead to a singular behaviour of the partonic cross section that becomes stronger with 
every order in The consequence is that the corresponding hadronic cross 

section does not converge as successively higher orders in 1/Mt an included. Our 
solution of this problem will be described below. 

• For the gg channel, the coefficient of the 1/x term at i = 2 turns out to vanish. 
The strategy of Ref. [35] for analysing the 1/Mt corrections at NLO was therefore to 
discard the terms with z > 2 in the gg channel, and to replace the 1/Mt expansion for 
the other channels by the exact result (which is much easier to obtain than for the 
gg channel). For the NNLO case, we do not have that option, because the calculation 
of the full mass dependence is currently out of reach. As mentioned before, we will 
discuss the problem of low x (large s) in more detail in Section 5. 

Let us now turn to the NNLO results. For the soft+virtual terms, we obtain: 



(2) 11399 19, 133, 165, 9, / 1189 2, 5, 5,' 
^ ' + -Trim + — C2 - — Cs - 7:C4 + ni --— + -Im - -(2 + 7C3 



" 144 8 2 " 4 ' 8 V 144 

(2) 47437199 /89 7 , \ 1909181 ^ 883 
+ C2 77 + 77 + ..... Cs + ttttt 



,(2) 



1244160 " \45 45 J 55296 " 1080 

/ 14563 281 
^ V 48600 ~ 2880 ^* ~ 90 , 

998645169149 , / 9677 857 , \ 267179777 , 4039 , 



117050572800 ^ V 37800 37800 J 35389440 ^ 51840 
/ 4565713 857 193927 
^' V 285768000 ~ 75600 ~ 21772800 ^* 
(2) _ 1712964005499545249 / 646571 17881 ^^^\ 5756378217151 
~ 39328992460800000 ^ V 15876000 ^ 4536000 j + 158544691200 
88077779 / 8432587511 17881 111726613 

^ 7620480000 ^* ^ V 4800902400000 ~ 9072000 '"^ ~ 91445760000 ^* 

^ 3 45 ' ^'^ 9450 ' ^'^ 15876000 ' 

6^^) = -33 + 2 , b'-^J =0 V i > 2 , 

45 = 72, 45 = Vi>2, 



(18) 



where Cn = C('^) is Riemann's zeta function with values 



2 4 

C2 = ^ = 1.64493... , Cs = 1.20206... , C4 = ^ = 1.08232 ... , (19) 
6 90 
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and 



(20) 



The various hard contributions are evaluated as expansions around x = 1. We quote them 
through 0(1 — x) for the sake of brevity. Terms through order (1 — x)^^, for arbitrary 
values of fip and /iR, are available upon request from the authors. 

For the gg channel, we find 



1453 
12 " 

+ ni 

+ {l-x) 

2385 



147C2 -35IC3 



77 1 01 

-+AC2 + -^L{x)-4L'{x) 



3437 1017 
— ^ + 



^ , , /2379 



1053 
T 



270 C2 



h 



(2) 



h 



h 



.9/ N N / 395 22 ^ 45 , 22 , 

L^(x) + 216 L^(x) + ni C2 L(x) + — L^(x 

4 ^ ^ ^ ^ \ 24 3 2 ^ ^ 3 ^ ' 



68 272 



99,2 45 45 



+ (1 - x) 



6661 9 ^ 172 , 

C2 H L(x) 

1200 80 ^ 15 ^ ^ 



80 ^ ' 



27 
40 



+ ••• , 



Al 1 r/ N 

+ re; L[x) 

' V80 20 ^ ^ 

3553 3553 



+ 



(2) 

^^•■^ 18900 4725 



Lix) + {l-x) 



1837333 21 
10584000 ^ 3200 '"^ 



189 ,2/ N 



471 



3200 
917641 



11200 
917641 



Im + ni 



659 7 , , 



49927 
50400 

+ 



L(x) 



(2) 

^^•^ 31752000 7938000 



L{x) + (1 - x) 



67200 2400 



13461173 1697 ^ 5644979 , 



635040000 896000 



+ 



15273 
896000 



L^x) + 



15731 
1344000 



Im + ni 



23347 1697 ^ 
8064000 ^ 2016000 



42336000 

+ 



(21) 



Here and in the following equations, the ellipse indicate higher orders in (1 — x). The qg 
(and qg) channel reads 
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,(2) 



^99,2 



(2) 



''99,e 



11 29 



311 



85 



367 



\81 3 ^ ^ 18 ^ 7 ^ ^ V 18 9 ^ 



+ (l-x) 



959 „ ^ 433 ^ , , 

h8C2H L(x) 

18 ^ 9 ^ ^ 



68 136 , , , 
405 + 405 ^(^) + (^ 



62737 
24300 



33 

' Y 

539 



1441 , / 44 11 ^ , , 

3240 V405 270 ^ ^ 



3553 



3553 



''5'^ 170100 85050 



L{x) + {1 



1080 

+ ••• , 

66227323 



C2 + 



4 
9 

4367 



L^{x) + -niL{x) 



+ 

187 



_ 7157 
~ 302400 
917641 



L\x) 



+ 



37481 
2721600 
917641 



Im + ni 



285768000 
421 



3240 360 



2947 ^ 26401 , 
^2 + 7;^;e^L{x) 



(22) 



129600 
421 



285768000 142884000 



L{x) + {1 



85050 226800 
1018432391 



L{x) 



388800 



+ ••• , 



39011 ^ 2460599 ^, , 

C2 + oo^.^oooo -^(^) 



34292160000 
94741 ^ 96389 



15552000 

/ 5573 
+ ni[ 



326592000 
5573 



36288000 



L'{x) 



65318400 



10206000 27216000 



Lix 



The remaining channels only start at higher orders in (1 — x): 



'"99,0 


— "'qqfi 


- %q'fi 




%q,i 


- %q',i 



(1-3 
0((1 



20 

xf) 



V i > 2 



16 ^ , , 16 \ 

— L(x) H L^(x) 

9 ^ ^ 9 ^ ^ 



+ 



(23) 



Again, we have included the known leading terms in 1/Mf for the sake of completeness. 



5 Small-x behavior 



As pointed out above, the 1/Mt expansion cannot give the proper result in the low-x 
(large-s) region, which is why the (1 — x)-expansion in this approach is no worse than 
the full X dependence. Fortunately, in Ref. [49], it was derived that the leading small-x 
behavior of the partonic cross section is given by (/xp = fJ-R = Mu) 

aW(x) =3(JoC(i) +0(x), <Tg)(x) = -9aoC(2)lnx + c + 0(x), (24) 

where the coefficients C^^'^ and C^^^ are available in Ref. [49] in the form of a numerical table 
for various values of Mt/Mn out of which we constructed simple interpolating functions. 
The constant c was undetermined. Note that Eq. (24) is the actual limit of s — > oo, i.e., it 
is not derived in the heavy-top limit. 
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We may use this additional information to improve our result in the following way: 



where agg denotes the expansion of the partonic cross section around x = 1 through 
0{{l — x)^). Note that ag^J{x) and afg{x) have the correct behavior for x — > and x — > 1 
up to the orders considered. In this way, we arrive at smooth functions that approximate 
the full partonic cross section over the full x-range. In order to illustrate the quality of 
this method, Fig. 3 shows (^^^^(x) together with the soft expansion CF^gg'^ {x) for iV = 8, as 

well as cr^gg,oo, i-e., the full x-dependence of the heavy-top result. All expressions include 
top mass corrections through 0{1/Mi); the curves are normalized by uq. 

In order to be able to directly compare Fig. 3 with Ref. [50] (which updates the numerics of 
Ref. [49]), the scales in Fig. 3 (a) are chosen identical to those of Figs. 1 and 2 in Ref. [50]. 
For the same reason, we set Mt = 170.9 GeV at this point, in contrast to the actual 
numerical section of this paper where the current world average for Mt is used. And 
finally, the curve for cr^gg,oo in the heavy-top limit is included (i.e., only the l/M^ terms), 
as in Ref. [50]. In case one is misled by the apparently large effect, Fig. 3 (b) shows the 
same curves on a linear scale. In fact, as will become obvious shortly, the effect of this 
matching on the hadronic cross section is rather small compared to using just the pure 
soft expansion ag^J'^{x) (cf. Fig. 5 below). 

Similarly, Fig. 4 shows the NNLO partonic cross section (gluon-gluon channel) as con- 
structed from Eq. (25), again with choice of scales and set of parameters as in Ref. [50]. 
The agreement to the result obtained in Ref. [50] is good, with a small difference on the 
unknown constant c in Eq. (25). The effect of the matching on the hadronic cross section 
is completely negligible when compared to the pure soft expansion a^gg'^ {x) (cf. Fig. 7 be- 
low). This is of course due to a suppression of the large-x region by the parton densities. 

6 Numerical results 

Unless stated otherwise, in all our numerical analyses in this paper, we use Mj = 173.1 GeV 
and set /ip = = -^H- The latter restriction shall be sufficient for this first study of the 
top mass effects; more detailed phenomenological studies, including scale variations, are 
left for a future publication. 

6.1 Next-to-leading order 

The natural extension of the heavy top limit of Eq. (8) would be to use Eq. (2) with the full 
top mass dependence in (Tq and the l/M^ expansion of A^^ as given in Section 4, and to 
match the result to the large-s region as defined in Section 5. However, in order to strictly 
test the heavy-top limit, we prefer to apply a consistent 1/Mt expansion to the partonic 

^The upper equation essentially corresponds to the matching procedure suggested in Ref. [49]. 





4y'^(x) + (l-x)^+i [3.oC«-a«'^(0) 

r ^1 1 



(25) 



fc=i 



12 



A A 




(a) (b) 

Figure 3: The NLO contribution to the partonic cross section function 
of X = My/s on a logarithmic (a) and a hnear scale (b). The expansion in (1 — x) 
converges up to the threshold x = Ml/{4M^) « 0.14. This expansion, through 

(1 — x)^, is displayed as the dashed curve, <7gg^'^, N = 8. The full x-dependence of 
the l/Mf expansion is shown as the dotted curves (lower: leading term in M}i/Mt, 
upper: including terms of order 1/M^). The solid line is the combination of the 
(1 — x)-expansion with the leading behavior at x ^ 0, cf. Eq. (24). 




(a) (b) 

Figure 4: Same as Fig. 3, but at NNLO. The dotted line is the full x-dependence 
of the 1/Mj° result. 
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cross section, without factoring out the LO mass dependence. Once the convergence of 
this 1/Mf expansion and its consistency with the heavy-top hmit of Eq. (8) is shown, one 
could try to derive an "improved heavy-top hmit" by keeping the fuU mass dependence 
in (Tq. However, as we will see, the improvement achieved in this way is well below any 
expected experimental accuracy. 

At NLO, we therefore define 

^T(^r) = ^0 s^ghA^ - ^) + V (^") ' (26) 

where d'l^^{Mp) is the NLO contribution to the partonic cross section evaluated as an 

expansion through 0{{M}i/AIt)'"'). It is obtained by expanding gqA^^^ with ao and A^^^ 
from Eqs. (3), (9), (12)-(16) in terms of 1/Mf, and applying the matching procedure 
of Eq. (25). The corresponding hadronic quantity derived from Eq. (26) is denoted by 
(T^^^(A/"), as usual. Note that it also depends on the depth of the expansion in (1 — x). 

First, we look at the convergence of the 1/Mj expansion of the gg channel alone, whose low- 
X behaviour is implemented as described in Section 5. Fig. 5 shows the ratio of ct^^^(M"), 
keeping various orders in 1 /Mt , to the fully mass dependent result C7™*^^u (dashed: l/M^", 
n = 0, . . . , 8; solid: l/M^^). The dotted line corresponds to the pure soft expansion result 
a^g, without matching to the low-x behaviour. The convergence towards the exact result 
is excellent, both for the LHC (a) and the Tevatron (b). The slightly better behaviour for 
the Tevatron is due to the smaller high-energy region. The effect of the matching from 
Section 5 is rather small. 

Unfortunately, the analytic low-x behaviour is currently not known for the subleading 
channels {qg, qq, and also the NNLO channels qq, qq'). However, their numerical contribu- 
tion at NLO is very small (for "reasonable" values of ^-p ~ Mh): for Higgs masses between 
100 and 300 GeV, the qg contribution is at the few percent level (< 4%) and the qq at the 
permille level. 

One could therefore just proceed by ignoring the 1/Mt-corrections at NNLO for all but the 
gg channel. At NLO, however, we observe that we can nevertheless improve the prediction 
by including the l/Mt and (1 — x) expansion of the subleading terms up to a certain 
depth beyond which the convergence properties of the series for the individual channels 
deteriorate, as is typical for an asymptotic series. Applying this criterion at NLO allows 
us to include the terms through order for the qg channel, but only the first two 

terms for the qq channel. The convergence of the expansion in (1 — x) is excellent (see 
below). 

This "optimal" result will be denoted by a^^^{Mt). It is shown in Fig. 6, divided by the 
exact mass dependence, both for a 14TeV proton-proton collider (a), and for a 1.96 TeV 
proton-anti-proton collider (b). One observes a nice convergence for the soft expansion 
towards a result that reproduces the full mass dependence at the sub-percent level. For 
comparison, the figure also contains the ratio of the heavy-top limit to the exact mass 
dependence (dotted line) which deviates from one at the percent level, as pointed out 
before. 
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Figure 5: Ratio of the gg induced component of the NLO hadronic cross section as 
obtained from Eq. (25) to the value obtained from HIGLU [51], when keeping suc- 
cessively higher orders in 1/Mt (decreasing dash-length corresponds to increasing 
order); the dotted line is the result obtained from the pure soft expansion <7gg 
through order without the matching of Eq. (25). 
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Figure 6: Solid: the final result for the NLO cross section, divided by the full top 
mass dependent result as obtained from HIGLU [51]. The dashed lines correspond 
to various orders in the (1 — x) expansion and nicely demonstrate the quality of 
the convergence. The dotted lines represent the result of the heavy-top limit at 
NLO (Eq. (8)). 



15 



6.2 Next-to-next-to-leading order 



In analogy to Eq. (26), we define 



i.NNLO,' ji/f" 



TT 



af3,oo 



as 

+ ( — 

vr 



CJ 



(27) 



where (T^^(M") is obtained by expanding a^A^^^^, with (Tq and A^^^ from Eqs. (3), (9), 
(18)-(23), in terms of 1/Mt, up to power n, and applying the matching procedure of 
Eq. (25). We already know that the first term in Eq. (27) is an excellent approximation to 
the full top mass dependent result, and therefore Eq. (27) provides a suitable quantity to 
compare the heavy-top result of Eq. (8) to the 1/Mt expansion. 

Again, we first look at the convergence of the 1/Mt expansion of the gg channel alone, 
whose low-x behaviour is implemented as described in Section 5. Fig. 7 shows the ratio 
of cr^^^*^(M"), keeping various orders in 1/Mt, to the heavy-top result o'^^^ of Eq. (8), 
which we recall includes the exact LO mass dependence (dashed: 1/M", n = 0,2,4; 
solid: 1/AIt). The dotted line corresponds to the pure soft expansion result a^g, without 
matching to the low-x behaviour. We observe very good convergence towards the heavy- 
top result, assuring us of the high quality of the latter. 

We then apply the same criteria as at NLO in order to obtain the "optimal" result 
(T^^^^(Mj) for the Mt dependent NNLO terms. They allow us to keep all four avail- 
able terms for the gg and the qg curve, i.e., through order 1/M^, and again only the first 
two terms for the qq, qq, and qq' initiated sub-processes. The result, divided by the one 
obtained from Eq. (8), is shown in Fig. 8, for various orders in the soft expansion which 
again converges excellently. 

The final result is within 0.5% of the well-known heavy-top limit, clearly justifying its 
successful and extensive application in the literature. As outlined above, one can now try 
to improve the NNLO prediction by factoring out the LO top mass dependence. However, 
it is obvious that this will not alter the result by more than 0.5% which is way below any 
expected experimental accuracy. 



(2) 



,(2) 



7 Conclusions and Outlook 

Top mass effects on the NNLO Higgs cross section in gluon fusion have been calculated. 
For the dominant gg channel, the result for s < AMf has been matched to the limiting 
behaviour at x — > as obtained in Ref. [49] . We have demonstrated the reliability of our 
result through the excellent convergence of the soft expansion. The main result is the 
confirmation of the remarkable quality of the heavy-top limit as defined in Eq. (8) which 
agrees with the 1/Mt expansion to better than 0.5% in the phenomenologically interesting 
mass range 100 GeV < Mh < 300 GeV, both at the LHC as well as at the Tevatron. This 
is an extremely comforting result because it validates the numerous higher order analyses 
that have been carried out up to now, in preparation for the LHC experiments [52, 53], 
and in particular for the ongoing Higgs searches at the Tevatron [54]. 

It remains to be seen to what extent this finding carries over to less inclusive quantities like 
distributions or phase space cuts. This requires a substantial extension of our approach 
and is left for future studies. 
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Figure 7: Ratio of the gg induced component of the NNLO hadronic cross section 
as obtained from Eq. (25) to the heavy-top result of Eq. (8), (decreasing dash- 
length corresponds to increasing order in l/M^); the dotted line is the result 
obtained from the pure soft expansion cr^gg'^ through order l/Mt" without the 
matching of Eq. (25). 
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Figure 8: NNLO version of Fig. 6, but now the curves are normalized to the 
heavy-top limit of Eq. (8). 
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